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Abstract 

We determine the most general form of the interaction between the gravitational field 
and an arbitrary Yang-Mills system of fields (massless and massive). We work in the 
perturbative quantum framework of the causal approach (of Epstein and Glaser) and use 
a cohomological definition of gauge invariance for both gauge fields. We also consider the 
case of massive gravity. We discuss the question whether gravity couples to the unphysical 
degrees of freedom in the Yang-Mills fields. 
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1 Introduction 



The general framework of perturbation theory consists in the construction of the chronological 
products such that Bogoliubov axioms are verified [1], [3j, for every set of Wick monomials 
Wi{xi), . . . , Wn{xn) acting in some Fock space H generated by the free fields of the model one 
associates the operator Twi,...,Wni^ij • • • > ^n)', all these expressions are in fact distribution- valued 
operators called chronological products. Sometimes it is convenient to use another notation: 
T{Wi{xi), . . . , Wn{xn))- The construction of the chronological products can be done recursively 
according to Epstein-Glaser prescription [Ij, ^ (which reduces the induction procedure to 
a distribution splitting of some distributions with causal support). These products are not 
uniquely defined but there are some natural limitation on the arbitrariness. If the arbitrariness 
does not grow with n we have a renormalizable theory. 

Gauge theories describe particles of higher spin. Usually such theories are not renormal- 
izable. However, one can save renormalizability using ghost fields. Such theories are defined 
in a Fock space H with indefinite metric, generated by physical and unphysical fields (called 
ghost fields). One selects the physical states assuming the existence of an operator Q called 
gauge charge which verifies = and such that the physical Hilbert space is by definition 
"Hphys = Ker{Q) / Im{Q). The space 7i is endowed with a grading (usually called ghost number) 
and by construction the gauge charge is raising the ghost number of a state. Moreover, the 
space of Wick monomials in Ti is also endowed with a grading which follows by assigning a 
ghost number to every one of the free fields generating Ti. The graded commutator dq of the 
gauge charge with any operator A of fixed ghost number 



is raising the ghost number by a unit. It means that dq is a co-chain operator in the space of 
Wick polynomials. From now on [•, •] denotes the graded commutator. 

In a gauge theory one assumes also that there exists a Wick polynomial of null ghost number 
T(x) called the interaction Lagrangian such that 



for some other Wick polynomials T'^. This relation means that the expression T leaves invariant 
the physical states, at least in the adiabatic limit. In all known models one finds out that there 
exists a chain of Wick polynomials T'^, T^^ , j'/j-^p^ _ _ _ such that: 



In all cases T'^'^, j'p-'^p^ qj-q completely antisymmetric in all indices; it follows that the chain 
of relation stops at the step 4 (if we work in four dimensions). We can also use a compact 
notation where J is a collection of indices / = [z^i, . . . , z/p] (p = 0, 1, . . . , ) and the brackets 
emphasize the complete antisymmetry in these indices. All these polynomials have the same 
canonical dimension 



dQA=[Q,A] 




[Q,r] = ^9^T^ 




[Q, T] = zd,T>^, [Q, T>^] = ^9,T'^^ [Q, T'^l = id.T^^'"', . . . 
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and because the ghost number of T = T*^ is supposed null, we then also have: 

gh{T') = (1.5) 
One can write compactly the relations (11. 3p as follows: 

cIqT' = I d^T"'. (1.6) 
For concrete models the equations (11.31) can stop earlier: for instance in the case of gravity 

Cohomology problems of the type (11.31) have been extensively studied in the more popular 
approach to quantum gauge theory based on functional methods (following from some path 
integration method) [2] (and refs. given there). In this setting the co-chain operator is non- 
linear and makes sense only for classical field theories. On the contrary, in the causal approach 
the co-chain operator is linear so the cohomology problem makes sense directly in the Hilbert 
space of the model. For technical reasons one needs however a classical field theory machinery 
to analyze the descent equations more easily. 

In some previous papers [7], [8] one of us has developed a machinery of analyzing systemat- 
ically the descent equations. Here we want to apply these methods for the interaction between 
gravitation and Yang-Mills fields. We consider the case of massless and massive gravity. The 
case of massless Yang-Mills fields was analyzed in [S] and here we provide the general case. We 
will use a convenient geometric setting for our problem presented in [7]. 

In the next Section we introduce the gauge structure on the various free quantum fields. In 
Section [3] we remind the reader the cohomology of the operator for Yang-Mills models and 
gravitation. Using this cohomology and the algebraic Poincare lemma we can solve the descent 
equations in Section H] and we determine the interaction between massless (or massive) gravity 
and Yang-Mills fields in the most general case. 

2 Free Fields of Spin 1 and 2 

We remind here some results and notations from [7] and [8j. The Hilbert space 7i. we use is of 
Fock type and it should describe particles of spin 1 and 2 with null or positive mass; we will 
denote by Q the vacuum state in H. The Pauli- Jordan distribution of mass m is denoted by Dm 
and Dm is its positive frequency part. The Minkowski metrics (with diagonal 1, —1, —1, —1) 
is denoted by r]^^. We will always mean by [■, ■] the graded commutator. The Hilbert space is 
generated by Bose and Fermi fields so it is useful to introduce a grading in Ti. as follows: every 
state which is generated by an even (odd) number of Fermi fields and an arbitrary number of 
Bose fields is even (resp. odd). 

2.1 Massless Vector Fields 

We consider a vector space Ti. of Fock type generated (in the sense of Botchers' theorem) by 
the vector field (with Bose statistics) and the scalar fields u,u (with Fermi statistics). The 
Fermi fields are usually called ghost fields and one can introduce in a natural way a ghost 
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number operator in the Fock space {u has ghost number 1 and u has ghost number —1). We 
suppose that all these (quantum) fields are of null mass. In this vector space we can define a 
masculine form < -, ■ > in the following way: the (non-zero) 2-point functions are by definition: 

< u{xi)u{x2)^ >= -i dI^\xi - X2) < ^, u{xi)u{x2)^ >= -i D^q^\xi - X2) (2.1) 

and the n-point functions are generated according to Wick theorem. To extend the masculine 
form to Ti we define the conjugation by 

Al = A^, = u, = -u. (2.2) 

Now we can introduce the operator Q according to the following formulas: 

[Q, A^] = I d^u, [Q, u] = 0, [Q, u] = -I d^A^ 

gi] = (2.3) 

where by [-, ■] we mean the graded commutator. One can prove that Q is well defined. Indeed, 
we have the causal commutation relations 

\A^,{xi),A^{x2)\=i ri^,^ Dq{xi- X2) ■ /, \u{xi)-,u{x2)\ = -% Dq{xx- X2) •/ (2.4) 

and the other commutators are null. The operator Q should leave invariant these relations, in 
particular 

[Q, [A^(xi), 'u(x2)]] + cyclic permutations = (2.5) 

which is true according to the previous definition. The usefullness of this construction follows 
from: 

Theorem 2.1 The operator Q verifies = 0. The factor space Ker{Q) / Ran{Q) is isomor- 
phic to the Fock space of particles of zero mass and helicity 1 (photons, gluons). 

Q is usually called gauge charge operator, the relations (12.31) define the gauge structure of 
the free fields. We notice that the operator Q raises the ghost number of a state (of fixed ghost 
number) by an unit. 

2.2 Massive Vector Fields 

We consider a vector space H of Fock type generated by the vector field A^, the scalar field 
$ (with Bose statistics) and the scalar fields u,u (with Fermi statistics). We suppose that all 
these (quantum) fields are of mass m > 0. In this vector space we can define a sesquilinear form 
< ■, ■ > in the following way: the (non-zero) 2-point functions are by definition: 

< Q, A^{xi)A^{x2)n >= i ri^^ D[^\xi - X2), < VL, $(xi)$(x2)fi >= -i D^^\xi - X2) 
< fi, u{xi)u{x2)VL >= -i DI^\xi - X2), < ^, u{xi)u{x2)il >= -i -D^^(xi - X2) (2.6) 
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and the n-point functions are generated according to Wick theorem. To extend the sesquihnear 
form to Ti we define the conjugation by 

Al = Af,, u} = u, = -u, = (2.7) 

Now we introduce Ti. the operator Q according to the following formulas: 

[Q,A^]=i d^u, [Q,u]=0, [Q, u] = (^^A^ + m $) [Q, $] = i m m, 

Qn = 0. (2.8) 

One can prove that Q is well defined. We then have: 

Theorem 2.2 The operator Q verifies = 0. The factor space Ker{Q) / Ran{Q) is isomor- 
phic to the Fock space of particles of mass m and spin 1 (massive photons, vector bosons). 

2.3 Gravitational Field 

We consider the vector space Ti of Fock type generated by the symmetric tensor field h^i, (with 
Bose statistics) and the (ghost) vector fields u^^u'^ (with Fermi statistics). We suppose that all 
these (quantum) fields are of null mass. In this vector space we can define a sesquihnear form 
< ■, ■ > in the following way: the (non-zero) 2-point functions are by definition: 

< VL, h^u{xi)hp„{x2)VL >= (r^^p r^j,^ + r]^p ri^„ - r/^j, rip„) d\^\xi - X2), 

< VL,Up{xi)u^{x2)^ >= i Vfj-iy -X2), 
< Q.,Up{xi)uu{x2)^ >= -i r]^y d\^\xi - X2) (2.9) 

and the n-point functions are generated according to Wick theorem. To extend the sesquihnear 
form to Ti we define the conjugation by 

= hpy, ul = Up, ul = -Ua. (2.10) 

Now we can introduce the operator Q according to the following formulas: 

[Q, hpu] = {dpU^ + duUp - r]pi,dpuP), [Q, Up] = 0, [Q, Up] = i d^hpy 

QVL = {) (2.11) 

where by [■, ■] we mean the graded commutator. One can prove that Q is well defined. Indeed, 
we have the causal commutation relations 

1 

[hpy{xi), hpa{x2)\ = -- {rjpp rj^„ + rj^^p r]p„ - rjpy rjp„) Dq{xi - X2) ■ /, 

[u{xi),u{x2)\= i r]py Dq{xi- X2) ■/ (2.12) 

and the other commutators are null. The operator Q should leave invariant these relations, in 
particular 

[Qi [hpui^i)^ u„{x2)]] + cychc permutations = (2-13) 

which is true according to the previous definition. The usefuUness of this construction follows 
from the following result [6]: 
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Theorem 2.3 The operator Q verifies = 0. The factor space Ker{Q) / Ran{Q) is isomor- 
phic to the Fock space of particles of zero mass and helicity 2 (gravitons). 

2.4 Massive Gravity 

We consider a vector space Ti. of Fock type generated (in the sense of Botchers' theorem) by 
the tensor field h^^, the vector field (with Bose statistics) and the vector fields u^,u^ (with 
Fermi statistics). We suppose that all these (quantum) fields are of mass m > 0. In this vector 
space we can define a sesquilinear form < •, • > in the following way: the (non-zero) 2-point 
functions are by definition: 

< fi, hf,^{xi)hp^{x2)^ >= -- (r/^p rjua + Vup V^t^ - Vi^u Vpa) -D^Ha^i - X2), 

< Q, u^{xi)uu{x2)Q >= i r]f,^ DI;![\xi - X2), 
< n,u^{xi)u^{x2)Vl >= -i 7]^^ Dl;l\xi - X2), 

< n,Vf,{xi)Vf,{x2)n >=i r]^^ ^^^xi -X2) (2.14) 

and the ra-point functions are generated according to Wick theorem. To extend the sesquilinear 
form to Ti we define the conjugation by 

= h^^, ul = Up, ul = -Ua, vl = v^. (2.15) 

Now we can introduce the operator Q according to the following formulas: 

[Q, h^y] = -- (dpUu + duU^ - rj^i^dpu''), 
[Q, M^] = 0, [Q, Up] = i {d^hpi, + mvp), 

[Q,Vp] = — — Up 

QVL = 0. (2.16) 
One can prove that Q is well defined. Indeed, we have the causal commutation relations 

[hpu{xi), hp^{x2)] = -- [ripp rj^„ + -q^p rjp^ - rjp^ rjp^) Dm{xi - X2) ■ I, 

[u{xi),u{x2)] = i r]p^ Dm{xi - X2) ■/ 
[Vp{Xi)Vp{x2)] = i Vpu Dmixi -X2) -I (2.17) 

and the other commutators are null. The operator Q should leave invariant these relations, in 
particular 

[Q, [hpu{xi), Ua{x2)]] + cyclic permutations = 0, 
[Q, [vp{xi),Ua{x2)]] + cyclic permutations = 0. (2.18) 

We have the result [9]: 

Theorem 2.4 The operator Q verifies = 0. The factor space Ker{Q) / Ran{Q) is isomor- 
phic to the Fock space of particles of mass m and spin 2 ( massive gravitons ) plus a spin 
particle of mass m. 
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2.5 The General Case 



The situations described above are susceptible to the following generalizations. First we consider 
the Yang-Mills case. We take a system of ri species of particles of null mass and helicity 1, that 
means we use in the first part of this Section ri triplets (A^, Ua, Ua), a G /i of massless fields; here 
Ii is a set of indices of cardinal ri. All the relations have to be modified by appending an index a 
to all these fields. In the massive case we have to consider r2 quadruples {A^, Ua, Ua, $a), a ^ h 
of fields of mass as in section 12.21 here I2 is a set of indices of cardinal r2. We want to 
include some arbitrary scalar fields with indices a G /s. Then we take / = Ji U /2 U /a a set of 
indices and for any index we take a quadruple Ua, Ua, $a), a G / of fields with the following 
conventions: (a) the first entry are vector fields and the last three ones are scalar fields; (b) 
the fields A^, are obeying Bose statistics and the fields Ua, Ua are obeying Fermi statistics; 
(c) For a G /i we impose $a = and we take the masses to be null nia = 0; (d) For a G /2 we 
take all the masses strictly positive: rria > 0; (e) For a E I3 we take A^,Ua,Ua to be null and 
the fields $a = 0^ of mass > 0. The fields Ua, a E hU I2 and a E I2 are called 
ghost fields and the fields 0^, a E I3 are called Higgs fields; (f) we include spinorial matter 
fields also i.e some set of Dirac fields with Fermi statistics: \E'a, A E h; (g) we consider that 
the Hilbert space is generated by all these fields applied on the vacuum and define in H the 
gauge charge operator Q according to the following formulas for all indices a E I : 

[Q, Ua] = -i (<9X + $a) [Q, $a] = i Ua, (2.19) 

[Q,^a]=0, (2.20) 

and 

Qn = 0. (2.21) 

If we want to include gravitons also then we extend the Fock space including the corresponding 
free fields /i^^, m^, u'^ (for massless gravity) or h^u, u'^,u'^, v'^ (for massive gravity) and we extend 
the definition of the gauge charge Q in a natural way using (12.111) or (12.161) . In this way the 
Fock space will describe a system of Yang-Mills particles together with (massless or massive) 
gravitons. 

3 The Cohomology of the Operator dq 

We know that the condition [(5,T] = i d^T^ means that the expression T leaves invariant the 
physical Hilbert space (at least in the adiabatic limit). 

Now we have the physical justification for solving the cohomology problem namely to deter- 
mine the cohomology of the operator dq = [Q, ■] induced by Q in the space of Wick polynomials. 
One can solve this problem in a quite general setting using the jet bundle formalism [7] and 
[8]. However for all practical purposes we need only a very particular case. Let us denote by 
•ptn,5 ^Y^Q set of Wick polynomials in the fields and their derivatives subject to the following 
restrictions: if T G P*"'^ then T is tri-linear in the fields (and their derivatives) and oj{T) < 5. 
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We are interested in the cohomology of the operator dq acting in P*''''^. We will denote by Zq' 

and B^q'^ the co-cycles and the co-boundaries of this operator, respectively. 

We introduce some notations. Basically we are looking for gauge-invariant variables. 
For the case of massless or massive spin 1 fields we define the field strength according to 

and observe that 

dgF^" = 0, 

dpF^u + d^F.p + d,Fp^ = (3.2) 

the last relation being called Bianchi identity (or homogeneous Maxwell equation by physicists). 
We denote by -Pj^-p the traceless part (in all indices) of the expression dpF^y (which also verifies 
the Bianchi identities). 

In the case of a massive vector field it is convenient to introduce another notation, namely: 

0^ = dp^-m (3.3) 

and we observe that 



V = 0. (3.4) 
We denote by cjyf!}, the traceless part of the expression 

dpd,^ - ^ m {d^A, + d,A^). (3.5) 



For the general Yang-Mills system we have to append an index a E IiU and we denote 
by P*"'^ C P*"'^ the space of polynomials in the variables Fg^p^^ Faf}u;p {o, G /i) and 
Fa^iu, F^%-p, (pa^i, (t>^a}v (o G h)- Then we have the following result [7]: 



Theorem 3.1 Letp G Z^q'^. Thenp is cohomologous to a polynomial of the form p = pi + dQP2 
where pi G p*"'^ and p2 G V^"'^. 

The expressions of the type dQP2 do appear because for massive Yang-Mills fields it is 
possible that the gauge charge operator does not raise the canonical dimension (for instance 
this is true for a monomial having a factor $a) 

In the case of the gravitational field it is also convenient to introduce some other notations: 
first ^ 

h = rj^^h^^ h^y = h^^- - r^^^ h (3.6) 

and the we define the Christoffel symbols according to: 

^fj.;up = dph^y -\- djjhpp d^hyp. (3-7) 

The expression 

Rpu;pa = dpTp.^^ -{p^ a) (3.8) 
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is called the Riemann tensor, we can easily prove 



Rfiu;p(7 Rvjj,;pcr Rfiu;(Tp Rp(T;iii/i 



dQRpu;pa — 0, 
Rpv\pa ~\~ Rpp]ua ~\~ R^a\up 0, 
d\Rf_iu-^pcr + dpR^y-fjX + da-Rpu;Xp — (3-9) 



the last two relations are called Bianchi identities. 

3(0) 
p,u;pa;X 



We denote by Rfl„„.x- the traceless part in all indices of d\R^pu;pa- and we also define 



Upu = M[H = ^ {dpU^ - dyU^). (3.10) 



Now we have the following result [8]: 



Theorem 3.2 Let p G Then p is cohomologous to a polynomial in u^, u^u and in 

p(0) d{0) 

In the case of massive gravity we also define the expressions 

^fiu = dpVy + duVp - iip^dpv'' - m h^,y 

(^ = r/^>^, (3.11) 

and observe that we also have 

dq^pu = 0. (3.12) 

By (^fS;p we denote the traceless part of dpif^i,. 

As in the Yang-Mills case we denote by P*"'^ C P*"'^ the space of polynomials in the 
variables Rpu-pa, R^^l-p^-x and Lpp^, ip, ^f}-p, dpip. Then we have [8]: 

Theorem 3.3 Letp G Z^q'^. Thenp is cohomologous to a polynomial of the form p = pi + dQP2 
where pi G "P*"'^ and p2 G P*"'^. 

We can obtain from the theorems above the description of the co-cycles in the general case 
of Yang- Mills fields interacting with massless or massive gravity using Kiinneth theorem [2]. 
This means that we must consider polynomials in invariants of both Yang-Mills and gravity 
type. 

We mention in closing that for polynomials from V^^^'^ Poincare lemma holds i.e. we have 
for any T^^ G P*"'^ 

5^T^ = =^ r^ = 9^T'^'' (3.13) 
where the expression T^^ G V^^^'^ is antisymmetric [7]. 
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4 The Interaction of Gravity with other Quantum Fields 



We here consider a system of massive and massless Yang-Mills fields {v^, Ua, Ua, ^a), a & I and 
we determine the coupling with the massless gravitational field hfj_u, u'', . By definition the 
ghost number is the sum of the ghost numbers of the YM and gravity sectors. As in [8] we 
consider that the interaction has null ghost number and its canonical dimension is bounded by 
5. We will get an expression of the form 

Ti rpYM I rpscalar i rp Fermi /a -t \ 

int — J int +-'int + int \^-^) 

with the scalar and Fermi contributions being the same as in [8]. We concentrate only on the 
Yang- Mills contribution and we have our main result. We use the following definitions. A Wick 
polynomial T is called a relatively co-cycle iff it verifies the relation dqT = d^T^] two Wick 
polynomials are relatively cohomologous iff the differ by an expression of the type (IqB + d^B^. 
As it was already seen in previous papers, one can reduce the problem of determining the 
relative cohomology groups to the cohomology of the operator dq using the descent procedure. 



Theorem 4.1 The expression T™ is relatively cohomologous to 

a,beli 

+ fab (4 V - h F,p^ Fi;'' + Au, d,u, Fr - 4 V € <Pb -^rriau^ 0^4-2) 

with the constants fab symmetric fab = fba and real, 
(a) The relation dqij^ = i d^t^^'^ is verified by: 

a,beli 

+ Y (""^ Fbp. + 4 F^-" F,,p - 2 <f)au K + ^mau, ^a <Pb) (4.3) 

a,bel2 

and we also have 

rfgC''^ = 0. (4.4) 



Proof: (i) By hypothesis we have 

rfQTint = ^ d.T.i^, (4.5) 
and the descent procedure based on Poincare lemma (see [7], [8] and [10]) leads to 

O^Q^int = ^ ^^^nt • 
"O-'int - ^ '^P^int 

dQTZ' = ^ d.TZ"'' 

dqTZr = (4.6) 
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and can choose the expressions T;^^ to be Lorentz covariant; we also have 

^7MT/J = |/|,^(T/J<5. (4.7) 

From the last relation and the theorems 13.11 and 13.21 we find that 

TZ'" = dQB^^^'^ + TZ'^^ (4.8) 

with TZ'q a polynomial in the invariants described in theorems 13.11 and 13.21 and we can choose 
the expressions B^^p" and T^'q completely antisymmetric. The generic form of T-Z^q can be 
easily obtained. If we substitute the expression of T-Z^" in the third relation (14.61) we find out 

dQiTZ' - ^ d^B^-n = ^ d.T-Z'o (4-9) 

so the expression in the right hand side must be a co-boundary and we immediately obtain 
T-Z'o = so: 

T-Z"'" = dgB"""" (4.10) 

and 

dQ{T.Z'-td^B'^^n=0- (4.11) 
We continue in the same way and obtain: 

= dgB^'^P + i d^B^^P" (4.12) 

and 

T-Z = dqB^"^ + I dpB'^'P. (4.13) 
(ii) We substitute the expression of T^, in the first relation (14.61) and get: 

dQ{TZ-zd.Bn = (4.14) 

Now it is again time to use theorems 13.11 and 13.21 to obtain 

TZ = dQB^ + zd^B>^'' + TZ^, (4.15) 

where q is a polynomial in the invariants appearing in the theorems 13.11 and 13. 2[ We get 
from the first relation (14.61) 

dgiT,^, - ^ d.B'^) = ^ d,TZ,o (4.16) 

so the right hand side must be a co-boundary. At this stage of the computation some non-trivial 
co-cycles do appear in o We consider only the Yang-Mills sector; then we have: 

^YM,/i _ ^(1) u T^pa 77 , fC^) P TPP.V 77 

-'int,0 ~Jab " bpa J " ^bup 

€^^^au4>l^9^iu,ra4>l^--- (4.17) 

where by ■ ■ ■ we mean other terms; we can impose the symmetry conditions 

/i;^ 9^^ = ^'^^- (4.18) 
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If one computes the divergence d^T^^^^^^ and imposes the condition that it is a co-boundary, 
then one gets: 

^(2) _ 4 .(2) (2) _ _2 (1) 

Jab ~ ^ Jab ^ iJab ~ ^ab ^ 

2 /i? + = (Va, beh), /i;^ = (Va G h, b G h) (4.19) 

and also the terms ■ ■ ■ from the expression of T^^q^ are null. The first equality above comes 
from terms without a mass factor and the second one from terms with a mass factor. It follows 
that one can take 

with tj^''^ the expression from the statement of the theorem. Because we have by direct 
computation dqtjj^ = i d^t(^'^ we get 

dQ{T,^,-t(^-id,B^^) = Q (4.21) 

so known results lead to 

^int = + dQB + i d^B'^ + Tint,o (4.22) 

where Tint,o is a polynomial in the invariants appearing in the theorems 13.11 and 13. 2[ But there 
are no such expression i.e. Tint,o = and we have 

T,^t = tlf + dQB + td^B^' (4.23) 

which is the final result. 

We point out that in all these computations one should consider all polynomials in the 
invariants, even or odd with respect to parity invariance. Indeed, because the Yang-Mills 
interaction is not parity invariant, there are no reasons to suppose that the interaction with 
gravity is parity invariant. Fortunately, the odd sectors do not produce non-trivial obstructions 
to the descent procedure. ■ 

Let us note that the result of the theorem stays true if we replace massless gravity by 
massive gravity. In the proof we will have to use now theorems 13.11 and 13.31 

We close by mentioning that the interaction between massive Yang-Mills fields and gravity 
(i.e. the second line in the formula (14. 2 p can be put into a simpler form. For simplicity we first 
consider one massive vector field i.e. I/2I = 1. In this case we can skip the index a = 1 and we 
define the physical part of the vector field according to the formula pi] : 

Af^' = A^ + \d^ A,. (4.24) 

This field has the following properties: 

dqAf^' = 0, d'^Af'^' = 0. (4.25) 
Then one can prove by some computations the following formula: 

h^, ra <Pl + ma Ua 0^ - d,Ua = h^" Af^' Af^' + dqB + 9^5^ (4.26) 
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i.e we can express the interaction between the massive vector field and gravity in terms of 
the physical part of A^ in a standard form: 

tint = T'^" (4.27) 
where we have defined h and h^i^ in formula (13. 6p and 

^ F,p F/ - - r/^, F,. F"'^ - Af^^ Af^^ + — r/^, A^'^^^'P Af^^ (4.28) 

is the energy-momentum tensor^ one can prove directly that it is conserved: 

T^" = 0. (4.29) 

In this form gravity couples to the physical degrees of freedom only, because there is no 
coupling to the ghost and A^^^ contains the three transverse physical modes only [11]. However, 
the new expression of the interaction Lagrangian has canonical dimension 7 so for the purpose 
of perturbation theory it is better to work with the expression appearing in the theorem. 

Note that the ghost couplings like d^Ua Fj^^ do not contribute to the S'-matrix elements 
between physical states in arbitrary order. Nevertheless, these couplings are necessary for gauge 
invariance of the theory. On the other hand the pure Yang- Mills couplings like fabc Ua A'^ d^Uc 
do contribute to physical S-matrix elements in ghost-anti-ghost loops. 

If one considers that the tensor fab in the formula (14 ■2p is diagonal then the energy- 
momentum tensor is additive with respect to the vector fields and the formula (14.271) extends 
naturally. (This property of the tensor fah follows from gauge invariance in the second order of 
the perturbation theory; this will be done in another publication). 

5 Conclusions 

In classical general relativity one says that gravity couples to everything which carries energy 
and momentum. That means one must have an energy-momentum tensor which is conserved 
du = 0. In quantum theory such a tensor does not exist in general. For example, in massless 
Yang-Mills theory and even in quantum electro-dynamics the naive (free) energy-momentum 
tensor T^*^ = F'^^ F'^p — 1 rj^^'^ Fp„ F^" is not conserved because 9^ = cannot be assumed 
as an operator equation. In this situation it is hard to find the correct gravitational couplings 
by classical Lagrangian arguments. 

Fortunately, the requirement of gauge invariance dqT = i d^T^ is so strong that it deter- 
mines all couplings uniquely, if some natural additional properties are assumed. It is a surprise 
that the resulting gravitational couplings contain ghost fields because they do not "carry energy 
and momentum" . The paradox is resolved by observing that the ghost coupling terms do not 
contribute to S'-matrix elements between physical states. 
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